Critical Behavior of the Meissner Transition in the Lattice London Superconductor 
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We carry out Monte Carlo simulations of the three dimensional (3D) lattice London superconductor 
in zero applied magnetic field, making a detailed finite size scaling analysis of the Meissner transition. 
We find that the magnetic penetration length A, and the correlation length £, scale as A ~ £ ~ \t\~"i 
with v = 0.66 ± 0.03, consistent with ordinary 3D XY universality, i/xy ~ 2/3. Our results confirm 
the anomalous scaling dimension of magnetic field correlations at T c . 
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The discovery of the high temperature superconduc- 
tors has revived interest in the effects of fluctuations on 
the critical behavior of the superconducting transition. 
The Meissner transition of a bulk type II superconductor 
in zero applied magnetic field is the most basic case that 
can be considered. While it was originally thought that 
this transition was weakly first order Q , it is now gener- 
ally believed to be in the same universality class as the 
ordinary three dimensional (3D) XY model, except with 
the temperature scale inverted M ^). The argument is 
based |^-^| on two observations: (i) the ordinary 3D XY 
model can be mapped onto a system of sterically interact- 
ing loops with inverted temperature scale Xi oop oc 1/T, 
and (ii) the vortex loops of a fluctuating superconduc- 
tor interact with a screened Coulomb interaction, with 
screening length equal to the bare magnetic penetration 
length Ao. Assuming that the finite interaction length Ao 
of the vortex loops is not a relevant modification of steric 
(on site) interactions, the universality of the fluctuating 
Meissner transition and the 3D XY model follows 0. 
Early Monte Carlo (MC) simulations by Dasgupta and 
Halperin || of a lattice London superconductor model 
strongly supported this picture by making a qualitative 
comparison of the shape of specific heat peaks in the two 
models. 

Recently there has been renewed interest in, and 
controversy concerning, the nature of this transition. 
Kiometzis et al. Q, considering a dual formulation of 
the fluctuating Ginzburg-Landau (GL) model, have ar- 
gued that while the correlation length £ ~ di- 
verges with the same exponent v as the ordinary 3D XY 
model, vxy ~ 2/3, the renormalized magnetic penetra- 
tion length should diverge as A ~ with v' = 1/2 the 
mean field exponent. Herbut and Tesanovic || however, 
using an analysis of the GL model exact to all orders in 
perturbation, have argued that, due to the presence of 
an anomalous dimension, tja = 1, for fluctuations of the 
magnetic field, one must have v = v' . Using a one- loop 
renormalization group (RG) scheme, they further sug- 
gested the possibility that v < vxy [0- Bergerhoff et 
al. [l2| , using a non-perturbative RG flow analysis of the 
GL model, similarly find v < vxy. Herbut plj however 



has argued that for the lattice London limit of the GL 
model, v = v' = vxy- 

To investigate this controversy we present here the re- 
sults of new MC simulations of the 3D isotropic lattice 
London superconductor (LLS) in zero external magnetic 
field. Carrying out the first detailed finite size scaling 
analysis of this model, we find results consistent with 
a single diverging length scale, hence £ ~ A. We find 
v ~ i/xY consistent with the universality of the ordinary 
3D XY model. We find clear evidence for the anoma- 
lous dimension of magnetic field correlations predicted 
by Herbut and Tesanovic ||. 

The Hamiltonian of our model 101 is 
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The sum is over all bonds of a 3D simple cubic lattice 
of unit grid spacing. 9i is the phase angle of the super- 
conducting wavefunction on site i, tpi = e l8i , where the 
amplitude of ipi has been taken constant (the London ap- 
proximation) . is the discretized vector potential on 
the bond at site i in direction (1 = x,y, z, and if /x, v, a is 
a cyclic permutation of x, y, z, then 



[D x A],-^ — Ai v + Aj_|_e,cr Aj 



A ia = 2nb Ul (2) 



is the counterclockwise circulation of the around the 
plaquette at site i with normal in direction /t. is the 
number of flux quanta 4> of total magnetic field through 
this plaquette. The coupling is J = 0§/167t 3 Aq, with Ao 
the bare magnetic penetration length, and U(ip) is the 
Villain function |H| 
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The first term in Eq. ([j]) is the kinetic energy of flow- 
ing supercurrents; the second term is the magnetic field 
energy. 

We focus here on the calculation of the magnetic field 
correlation function 
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where b^qv) = J^i &~ lqv ri b lti is the Fourier transform 
of the total magnetic field, and fi _L v. F{q) is just 
the wavevector dependent magnetic permeability, with 
\\m q -+ F{q) = dB/dH @. Our goal is to show that 
the singular part of F(q) is consistent with the scaling 
Ansatz 



F(t, Q, L) = l- x F{U x l v , Qt,L/l) 



(4) 



where t = T — T c , Q = 2 sin(g/2), L is the system length, 
and £ is an arbitrary length rescaling factor. Note that we 
choose Q rather than q as our scaling variable, since the 
vortex line interaction that arises from the Hamiltonian 
(|lj) is a function of q^ only through the combinations Q M 
fuj . Since Q — > q as q — > 0, this does not affect the long 
length scaling; our hope is that by using Q we may suc- 
ceed to slightly extend the scaling region to shorter length 
scales. Verification of the scaling Eq. (Q) will demonstrate 
that there is only a single diverging length scale in the 
model, that describes both the critical behavior of global 
thermodynamic variables, as well as the spatial varia- 
tion of magnetic field fluctuations. Since the former is 
determined by the correlation length £, while the later 
is determined by the magnetic penetration length A, we 
conclude that £ ~ A ~ 

We carry out standard Metropolis MC on the Hamilto- 
nian (|]) for cubic lattices of lengths L = 8 to 32, using pe- 
riodic boundary conditions. We use the particular value 
of Ao = 0.3 (in units of the grid spacing). Temperatures 
will be measured in units of J. In one MC "pass" we first 
update Ai X , Ai y and Ai Z at each site i, going sequentially 
through the entire lattice, then follow this by a sequen- 
tial update of the B. L . The Ai^ are allowed to fluctuate 
without constraint. For our largest system size, L = 32, 
we use at each temperature typically 32,000 passes to 
equilibrate, followed by 1.7 x 10 7 passes for computing 
averages. 

First we consider the scaling behavior of the magnetic 
permeability. Evaluating Eq. (Q) at the smallest wavevec- 
tor in our system, q min — 2ir/L, using Q m i n w q m i n , and 
choosing the rescale factor I = L, we arrive at 



LF(t, Qmin, L) = F(tL}t u ,2n, 1 



(5) 



Exactly at T c (i.e. t = 0), LF(q m i n ) should thus be a con- 
stant independent of L. In Fig.^Ja we plot our data for 
LF(q m i n ) vs. T, for L — 8 — 32. To a very good accuracy, 
the curves for different L do indeed intersect at a single 
point, T c w 0.8. To further verify the scaling relation 
Eq. (||), we fit our data for LF(q m i n ) near T c to a low or- 
der polynomial expansion in (T—T C )L 1 / 1/ . We determine 
the values of T c = 0.8000 ± 0.0002 and v = 0.66 ± 0.03 
from a 5th order polynomial fit, restricting data to the 
ranges |t| < t max = 0.006 and L > L min = 12. Increasing 
either the order of the polynomial, L m ; n , or decreasing 
i max resulted in no change in these fitted values, within 
the estimated statistical error. In Fig. |l|b we use these 
fitted parameters to plot LF(g m i„) vs. tl}l v ', for data in 



the range \t\ < 0.01. The resulting data collapse is very 
good. Our value of v is thus completely consistent with 
z^xy ~ 2/3. Note that by taking q — q m - m , L — > oo, and 
£ = t~ v in Eq. (|J), our results imply that the magnetic 
permeability vanishes, as T — + T+, as dB/dH ~ t v . 
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FIG. 1. a) LF(q) vs. T for system sizes L = 8 - 32. The 
common intersection of all curves locates T c . b) Scaling col- 
lapse of LF(q) vs. (T — T C )L 1 / 1 ' using fitted values of T c and 
v. 

We now consider the q dependence of F(q). In Fig.|2| 
we plot F(q) vs. Q for L — 8 — 32, exactly at T c and for 
one representative temperature above and below T c . We 
see, as expected, that for T > T c , F(q) approaches a con- 
stant as Q — > 0, while for T < T C) F(q) vanishes as Q 2 . 
Exactly at T c however, F(q) appears to vanish linearly 
as Q. This is a clear suggestion of the anomalous dimen- 
sion of magnetic field correlations predicted by Herbut 
and Tesanovic, according to which at T C) F(q) ~ q qA 
with r/A = 4 — D in D dimensions || . It is interesting to 
note that, while there is a considerable finite size effect 
for T > T c , finite size effects at a fixed value of Q appear 
negligible for all T <T C . 
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FIG. 2. -FX?) vs. Q for system sizes L = 8 — 32, at 
T = 0.76 < T c , T c = 0.80 and T = 0.82 > T c . Note the 
virtual absence of finite size effects for T <T C . 

To further verify the anomalous scaling dimension of 
magnetic field correlations, we can apply Eq. (Q) at t = 0, 
taking as the rescaling factor i = L, to get 



LF(0,Q,L) = F{0,QL,1) 



(6) 
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In Fig.g we plot LF(q) exactly at T c vs. LQ, for L = 
8 — 32 and q < n/2. We find a good collapse of the data 
to a single curve that vanishes linearly as LQ — > 0. 





1 




o 










- 


1 — 0.800 


O , '' 






q < tt/2 


o 






Py 










xL = 


8 - 




£ y ' 
p 


oL = 


12 " 






a L = 


16 " 






• L = 


24 " 




»'' 

i 


oL = 


32 ' 





10 20 


30 40 





12 
10 
8 
6 
4 
2 


U 111 zu 

LQ 

FIG. 3. Scaling collapse of LF(q) vs. LQ at T c , for 
L = 8 — 32, g < 7r/2. LF(q) vanishes linearly as LQ — > 0. 

Finally, in the thermodynamic limit L — > oo, we can 
use Eq. (g) with ^ = |*j — ^ = ^ to get 

F(t, Q, oo)/Q = (QO'^ifl, Q£, oo) , (7) 

where F± refers to distinct branches for T > T c and T < 
T c . Using the values of T c and v found in the fit of Fig. [l]b 
to determine £ = |T — T<\~ v , we plot in Fig.|] our data 
for F(q)/Q vs. £Q on a log-log scale. We use only data 
for which finite size effects appear to be small, and which 
are in the scaling region. We see an excellent collapse of 
the data. Fig. [| clearly demonstrates that there is only 
a single diverging length scale for the spatial variation of 
magnetic field correlations, and that this length scale is 
£. For the T > T c branch, we see that F(q)/Q diverges 
as as £Q — * 0, indicating that F(q) approaches a 

finite constant cx For the T < T c branch, we see 

that F(q)/Q vanishes as £Q as £Q — -> 0, indicating that 
F(q) vanishes as £Q 2 . However for both branches F(q)/Q 
approaches the same constant as £Q — > oo, indicating 
that F(q) vanishes linearly in Q exactly at T c . Fig. g thus 
gives another demonstration of the anomalous dimension 
of magnetic field scaling at T c . 

To get a better physical understanding of the effects of 
this anomalous dimension of magnetic field scaling, con- 
sider applying a small external magnetic field given by 
A ext . The London equation, describing the total screen- 
ing of the Meissner state, T < T c , gives for the induced 
supercurrent |15j 



7\2 



(8) 



where the vector potential of the total magnetic field is 
the sum of the applied and induced fields, (A^(q£>)} — 
tf?(qv) + {A™ d {qi>)), and 



\ 2 Q /a(q = 0) =n s (T)/n s (T = Q) 



(9) 
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FIG. 4. Log-log scaling collapse of F(q)/Q vs. £Q. The 
dashed lines at small £Q have slopes of ±1, to indicate the 
asymptotic behavior as £Q — > 0. 

The induced supercurrent is also related to (A™ A (qv)) 
by Ampere's Law, which for the gauge Q • A = can be 
written as 

(j^{qP)) = -J\lQ 2 {A i f{qi>)) . (10) 
Noting that F{q) is the magnetic permeability, we have 



F(q) = 



AT(qv) 



1 



Combining this with Eqs. (g) and (10) gives 

a{q)Q 2 



F(q) 



l + a(q)Q2 



(11) 



(12) 



is determined by the density of superconducting electrons 



Comparing with the results of Fig. g, we see that for finite 
£ at T < T c , we have lim^o a(q) ~ £. Eq. (^) thus im- 
plies that the superconducting electron density vanishes 
as n s ~ - \t\ v . 

The renormalized magnetic penetration length A is de- 
termined by the pole of F(q). If one could ignore the 
q dependence of ct(q), one would then conclude that 
A 2 = a(q = 0). From this follows A ~ and n s ~ A~ 2 . 
These are indeed the expectations from mean field the- 
ory Jig ] , as well as the "uncharged" superconductor rep- 
resented by the ordinary 3D XY model |ll| (given by the 
limit Ao — > oo). They also hold in the present model, at 
low temperatures. 

However, as T -> T~, Eq. @ and Fig.@ imply that 
lim g _>o a(q) ~ l/q. This is a consequence of the anoma- 
lous scaling dimension of the magnetic field. It is this 
singular dependence of a(q) on q that shifts the pole of 
F(q) so that A ~ £ rather than v 7 ?, in the "charged" 
superconductor critical region. In this critical region the 
London relation n s ~ A -2 no longer holds [p"7| . 

Note also that -F^q) determines the decay of the mag- 
netic field away from a test vortex line, hence it deter- 
mines the renormalized interaction between vortex lines. 
The q dependence of a(q) implies that in the critical re- 
gion, the interaction between two straight and parallel 
test vortex lines separated by distance r will change from 



3 



the In r of mean field (MF) theory at r < v 7 ? ~ Amf , to 
the faster decay of 1/r for < r < £ ~ A. 

Finally, we note that the anomalous scaling of F(q) 
also has some interesting consequences for the ordinary 
3D XY model. One can show that, within the mapping 
of the 3D XY model to a gas of sterically interacting 
loops, the helicity modulus of the XY model maps into a 
loop-loop correlation function. Identifying such loops as 
the vortex lines of the LLS, which as q — > (or Ao — > 0) 
become identical with magnetic flux, one concludes [ lq] 
that the wavevector dependent helicity modulus [ |T4jl9| 
T XY (g) of the ordinary 3D XY model should be the dual 
of F(q). We have carried out independent MC simula- 
tions of the ordinary 3D XY model, in the Villain ap- 
proximation, calculating T XY (q) for an ensemble with 
"fluctuating twist" boundary conditions (fbc) flj9| . We 
plot our results for T XY (q) vs. Q in Fig.|, for L = 8 - 32 
and T = 2.96 < T c , T c = 3.0 and T = 3.04 > T c . Note 
the striking similarity to Fig.||, only with the tempera- 
ture scale inverted. Finite size effects are negligible for 
T > T c . As Q — ► 0, T XY approaches a constant for 
T < T c , T XY - Q 2 for T > T c , and T XY - Q exactly at 
T c . Thus the anomalous scaling of F(q) at T c shows up 
as an anomalous scaling of T XY (q) at T c [E0[ . 
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FIG. 5. T XY (g) vs. Q for the ordinary 3D XY model. 
Data is for L = 8 - 32 at T = 2.96 < T c , T c = 3.0 and 
T = 3.04 > T c . Note the virtual absence of finite size effects 
for T>T C . 

To conclude, we have presented MC data that verifies 
the scaling Ansatz of Eq. (^). This Ansatz implies that 
there is only a single diverging length scale in the prob- 
lem, and that the magnetic penetration length scales as 
A ~ £ ~ |t|-". We find the value of v = 0.66 ± 0.03 con- 
sistent with vxy ~ 2/3 of the ordinary 3D XY model, 
and confirm the predicted anomalous scaling dimension 
of magnetic field correlations at T c . 
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